Breakdown of universal transport in correlated d-wave superconductors 
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The prediction and observation of low-temperature universal thermal conductivity in cuprates has 
served as a keystone of theoretical approaches to the superconducting state, but recent measurements 
on underdoped samples show strong violations of this apparently fundamental property of d-wave 
nodal quasiparticles. Here, we show that the breakdown of universality may be understood as the 
consequence of disorder-induced magnetic states in the presence of increasing antiferromagnetic 
correlations in the underdoped state, even as these same correlations protect the nodal low-energy 
density of states in agreement with recent scanning tunneling experiments. 
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Introduction. Thermal conductivity measurements at 
low temperatures T in the superconducting state have 
played an important role in strengthening the case for 
a d-wave BCS description of quasiparticles in optimally 
doped cuprate superconductors. They are bulk probes of 
the superconducting state, unlike angular resolved pho- 
toemission (ARPES) and scanning tunneling microscopy 
(STM), and can currently be performed at lower T 
than microwave experiments. One drawback is the need 
to separate phonon and electron contributions, but in 
the cuprate superconductors an asymptotic linear term, 
kq/T ^ const which can be attributed solely to quasipar- 
ticles, dominates at the lowest T. After theoretical pre- 
dictions of the universality of low T quasiparticle trans- 
port in nodal superconductors [1[, experimental confirma- 
tion was obtained in optimally doped materials [2]. Ac- 
cording to theory, which relies on the disorder-averaged 
self-consistent T-matrix approximation (SCTMA), the 
low-T thermal conductivity is 
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where ks is Boltzmann's constant, and wa, vf denote the 
nodal gap slope and Fermi velocity, respectively. This is 
a remarkable example of a transport coefficient which is 
unaffected-to leading order-by the addition of disorder. 
Furthermore, this result is insensitive to vertex correc- 
tions due to anisotropic impurity scattering and Fermi 
liquid effects (si . Observation of universal conductivity in 
optimally doped samples [J, |5| was a key step in establish- 
ing the existence of nodal quasiparticles and the validity 
of the SCTMA in this limit. Soon thereafter, however, 
measurements in underdoped samples exhibited values of 
the low-T thermal conductivity considerably below pre- 
dicted universal limits [6|. How and why the universal 
prediction ([T]) breaks down has never been explained. 

Since vf is generally considered to be well-known from 
ARPES [8], Eq.lJT]) has also been used to extract the 
gap slope for a number of cuprates at lower doping as 



well[6|, leading to the conclusion that va increases with 
underdoping[6|. This conclusion is in apparent contradic- 
tion to recent Raman[7| and ARPES experiments [8|, so 
it is even more important to examine physical effects out- 
side the framework of the SCTMA which could lead to a 
suppression of kq/T and thereby to a possible erroneous 
conclusion about the doping dependence of va ■ 

There are several effects known to lead to a suppres- 
sion of Kq/T. Localization effects were discussed in this 
context in Ref. 9] , and effects of bulk subdominant com- 
peting orders have been shown to suppress kq/T but 
do not immediately eliminate it, despite the removal 
of the d-wave nodes[10|. Here, we investigate the ef- 
fects on k{T) by local impurity-induced moments rele- 
vant e.g. to the underdoped regime of La2_a;Srj;Cu04 
(LSCO)[ll|, [H m. Such moments are formed in the 
presence of background Hubbard interactions in the host 
material [IJ]. 

Recently, the effect of strong correlations of this type 
on the density of states (DOS) of a disordered d-wave 
superconductor was investigated by Garg et al. [l5|, 
in an approach where the Gutzwiller approximation was 
used to approximately project out the doubly occupied 
states of a disordered t— J model. They found the rather 
remarkable result that the low-energy states were "pro- 
tected" by interactions, i.e. in the presence of the projec- 
tion the low-energy residual DOS which normally arises 
from disorder[l9| was strongly suppressed. Experimen- 
tally it is indeed observed that the low-energy DOS is 
surprisingly homogeneousfldllTI. Il8j . reflecting a robust- 
ness of the nodal quasi-particles to disorder. Garg et 
al. speculated that the projection mitigated effects of 
disorder generally at low energies. Here we show that 
in the density channel this hypothesis is correct: inter- 
actions screen and diminish the disorder potential. In 
the magnetic channel, however, scattering is enhanced, 
and properties which are sensitive to spin scattering like 
the universal thermal conductivity k{T) may be strongly 
renormalized even in the absence of Anderson localiza- 



tion. 

Model. The model we use to study disordered d-wave 
superconductors with magnetic correlations is 
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Here, c\^ creates an electron on site i with spin a, and 
tij = {t, i'} denote the two nearest neighbor hopping in- 
tegrals, y™P = ^ -^^ V™^P5ij is a nonmagnetic impurity 
potential resulting from a set of point-like scatterers of 
strength y*'"^^ ^i is the chemical potential adjusted to fix 
the doping a;, and Ay is the order parameter on the bond 
between sites i and j . The amplitude of Ay is set by the 
coupling constant g/t = 1.3. The Hubbard repulsion is 
treated by an unrestricted Hartree approximation, and 
U is taken to be homogeneous. Below we fix t' /t = —0.3, 
X = 0.1, and solve Eq.Q self-consistently by diagonal- 
izing the associated Bogoliubov-de GennesfBdG) equa- 
tions on systems of iV = 40 x 40 sites (20l.[21|. 

The model given by Eq.([2) has been used ex- 
tensively in the literature to study bulk competing 
phases, field-induced magnetization, as well as novel 
bound states at interfaces between antiferromagnets and 
superconductors [22]. It has also been used to study field- 
induced moment formation around nonmagnetic impuri- 
ties in correlated d-wave superconductors [201 l23l. |24J. In 
the case of many impurities, Eq.([2]) was recently used 
to model static disorder-induced antiferromagnetism as 
seen, e.g., by neutron scattering measurements |llL l2l|. 

Results. In the clean case (y*™P = 0), magnetic or- 
der induced by U will compete with the superconducting 
order and lead to a bulk magnetic state above a criti- 
cal value Uc2. For Uci < U < Uc2, a single point-like 
impurity can induce a localized S = 1/2 state with stag- 



gered magnetization [20, [23|, [2J, [25[. Here, Ud is the 
critical U necessary for local moment formation, and de- 
pends on band structure and impurity strength. When 
U < Uci there is no induced magnetization but the mag- 
netic correlations still suppress the charge modulations 
and shorten the superconducting healing length ^ near 
the impurity. 

The picture becomes considerably more complex for 
finite impurity concentrations. However, there is still a 
lower value U*^ below which no magnetism is induced, 
and an upper U*2 — Uc2 above which the ground state 
becomes that of a disordered magnet. As explained in 
Ref. [2l|, it easier for a dirty system to generate antifer- 
romagnetic islands {U*i < Ud), because of clustering of 
impurities and resulting charge redistributions. 

Density of states. We first focus on the DOS of the 
disordered d-wave superconductor in the presence of cor- 
relations and compare with the results of Ref. 15|. In 



particular, we would like to ascertain whether the onset 
of the local magnetic state has a qualitative effect on the 
total DOS. In Fig. [T^, we show the spatially averaged 
total DOS, Niuj) = J2r, N(ri,Lu)/N, where 

(3) 
for a series of disorder concentrations consisting of weak 
scatterers V^™P/t — 1.0 for the noninteracting case 
U = 0. We see the usual low-energy pile-up of impu- 
rity states inside the d-wave gap[26|. In order to obtain 
smooth DOS curves we have used a small artificial smear- 
ing factor 77/t — 0.025. Figure [TJj displays the results of 
the same study, but for the correlated d-wave supercon- 
ductor with U/t = 2.5. The remarkable result is that 
even for 25% disorder, the Hubbard correlations "pro- 
tect" the d-wave V-shaped D0S[l5|- The origin of this 
universal low-energy behavior is the suppressed charge 
modulations near the impurities as shown in Fig. [TJ;. 
Therefore, an impurity potential apparently disturbs a 
correlated superconductor much less than a conventional 
BCS state jl5.]. Interestingly, impurities which only mod- 
ulate the pair interaction will also perturb the electronic 
structure primarily near the antinodes, while protecting 
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FIG. 1: (Color online) (a) Spatially averaged DOS in the un- 
correlated d-wave superconductor ([/ — 0) for different con- 
centrations of disorder rii. One clearly sees the pile-up of low- 
energy impurity states, (b) Same as (a) but for the correlated 
d-wave superconductor with U = 2.5t. For a large range of 
impurity concentrations (0-25%) any sizable and qualitative 
change in the 'V'-shaped low-energy DOS is completely ab- 
sent, (c) A line cut through the center of the system showing 
the total charge density (top black lines) and the d-wave order 
parameter (lower red lines). The dashed (solid) lines display 
the results without (with) correlations, (d) DOS for 2.5% 
strong scatterers (1/""''/f = 10) as a function of U. The inset 
shows a zoom of the low-energy region uj £ [— 0.35t, 0.35t]. 



the low-energy DOS universality. This kind of disorder 
may exist in Bi2Sr2CaCu208+a; , where dopants seem to 
enhance the pairing interaction locallv[27ll28J. 

The results shown in Fig. [T^-c are for weak scatter- 
ers without induced magnetization, i.e. in the regime 
U < U*i. It is also interesting to study the DOS in the 
cluster spin-glass phase, U*i < U < U^, hi ghly rele- 
vant for example for underdoped LSCO Uj, ll2, |2l| . The 
tendency to create impurity-induced magnetic states in- 
creases with the impurity potential lJ| , so in Fig. [TJl we 
have compared the DOS with near-unitarity limit scat- 
terers with y""P/t =10 and a fixed impurity concen- 
tration of Tii = 2.5%. From Fig. [TJi it is evident that 
in the absence of correlations, the conventional plateau 
in the DOS predicted by the SCTMA is seen. As U is 
increased, this plateau disappears gradually, and there 
is no qualitative change in the DOS as U*i/t ~ 1.5 is 
crossed. Note at this relatively weak value of U, the pro- 
jection of doubly occupied sites has not been fully im- 
plemented, but the tendency is clearly the same. As U 
increases beyond Uc2/t ~ U*2/t ~ 2.8 we enter a regime 
where bulk magnetic order coexists with superconductiv- 
ity, and the clean DOS acquires an additional gap of the 
Mott-Hubbard type. However, as shown recently, even 
in this regime, disorder can result in low-energy DOS 



indistinguishable from a d-wave superconductor [29J. 

Thermal conductivity. For d-wave superconductors, 
a real-space BdG calculation of k(T) was used to 
study both localization phenomenajol, ll3| and the vor- 
tex state[30| . Here we follow the same approach but focus 
on the correlation effects and the role of disorder-induced 
magnetization using a realistic band structure for cuprate 
superconductors. For an inhomogeneous system we de- 
fine the total thermal conductivity k{T) as 

where Kxxi^i,T) = hx{vi)/{—V xT) is the ratio of the 
thermal current along the x axis at position r^ and the 
uniform temperature gradient applied along the x direc- 
tion. Within linear response we have 
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(5) 
where Qxx is the heat current-heat current correlation 
function. The quantity inside the parenthesis in Eq.([5]) 
reads in compact form 
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where the vector '0n(*) = [^tn(*): ^n(*)], and Vg denotes the discretized group velocity operator given by iig = Wfcj„T3 -t- 
'Sa'''iIO|- Here the kinetic velocity {v^in)ii = —Haixi — xi) with tu = i(i'/-\/2) for i,l being nearest (next- nearest) 
neighbors (in the flow direction), and similarly the gap velocity is given by {v^)ii = il^ii{xi — xi). In Eq.([B]) we have 
introduced the spatially independent thermal function F{En, Em) given by 
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where 5r^{u)) — 2rj/{Ld'^ + rj^) iQ\. The spin index (and the 
sum over spin) is implicit in Eqs.([n])-([7]). 

Numerical solution of Eqs.(in])-([7]) is computationally 
demanding, and we are restricted to systems of order 
40 X 40 sites. To study the low T regime, one needs to 
have a sufficient number of states within the d-wave gap. 
Therefore, for the discussion below we have increased the 
couphng constant to g/t = 2.8, giving Ay/f = 0.4 per 
link in the homogeneous case. This value of A/t gives 
a universal value kqo = 1/3 (wi^/'^A + i'A/f_F) — 0.88 
[ks — h = \) ior the band used here. We stress that 
the increase of the superconducting gap, which is a fac- 
tor 6-8 larger than cuprate materials, is merely a means 
to calculate k{T)/T at low T and is not important for 



the discussion below. However, the corresponding criti- 
cal [/'s will be larger as well (as discussed below). 

In Fig. [Ui (Fig. [5]d), we show the low-T thermal con- 
ductivity for the case of Ui = 2.5% (15%) scatterers with 
V^P/t = 10 (V'^^P/t = 3) averaged over 20 different ran- 
dom impurity configurations (which is enough for config- 
urational convergence). For U — the result in[2^ agrees 
well with those obtained previously for the dilute impu- 
rity limit: in the non-self-consistent (NSC) calculation 
(with homogeneous A=0.4), k(T)/T = kqo + aT^jSlI, 
whereas in the self-consistent (SO) case the spatial inho- 
mogeneity of A^j causes a reduction of a and leads to 
a more linear T dependence[9| of k/T. The case with 
ymip^-j- _ g Q ^^ [^j-, jg .^^ example of one of the weak- 
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FIG. 2: (Color online) (a-b) k(T)/T versus T as a function 
of increased U. (a) shows the case of 2.5% strong scatterers 
with V""''/t = 10.0, and (b) has 15.0% weaker impurities of 
yimp ^^ _ g g (^Q.(['j show the disorder- iuduced magnetization 
in real-space for the case in (a) with U/t — 3.2 (c) and U/t = 
3.6 (d). The inset in (d) shows the DOS (to G [-O.St, 0.3t]) for 
the clean system (dashed), and the disordered case from (a) 
with [7 = (red), U/t = 3.2 (green), and U/t = 3.6 (black). 



est disorder potentials we are able to study since fur- 
ther reduced T/*'"^ leads to mean-free-path larger that 
our system. As seen from Figs. [2^-b, an increase of 
the correlations U, initially increases k{T)/T towards 
the NSC result. The origin of this enhancement is the 
{/-suppression of the charge- and gap-modulations: the 
self-consistent mean fields approach the NSC result. As 
U increases further and magnetic scatterering centers are 
induced near the defects (Fig. [2fc-d), k{T)/T is contin- 
uously suppressed, and kq/T eventually vanishes in the 
bulk magnetic state. In BCS d-wave superconductors 
the universal ratio of kq/T arises from a cancellation 
of the increased disorder-induced scattering rate and a 
concomitant increase in the density of low-energy quasi- 
particles. In the present case, the universality is bro- 
ken because the disorder-induced moments increase the 
scattering rate whereas the resulting low-energy DOS re- 
mains unchanged. 

In the two cases shown in Fig. [2K-b, U*i/t ~ 1.9 
for V'"'P/t = 10.0, and U^^/t ~ 3.1 for y™P/t = 3.0, 
whereas Ucil^ = 3.9 in both cases. Density of states cor- 
responding to the magnetization plots in Fig. [2j;-d are 
shown in the inset of Fig. [2Ji. Thus, there exists a large 
spin-glass regime where the DOS remains universal in the 
sense that it maintains the characteristic 'V'-shape gap 
for a d-wave superconductor, but the thermal conduc- 
tivity is simultaneously suppressed due to the creation 
of additional effective magnetic scatterers. Therefore, 
impurity-driven local moment formation present in the 



underdoped regime may explain the doping dependence 
of k,q/T measured in LSCO 34|- Since the low-T thermal 
conductivity is strongly suppressed, a naive use of the 
"universal" clean d-wave result for kqo in Eq. ([T]) would 
lead to an erroneous estimate of the superconducting gap 
slope at the node. Thus the conclusion that ?;a is reduced 
with decreasing doping as found by recent Ramani?] and 
ARPES[3] measurements is not necessarily inconsistent 
with thermal conductivity measurements (6|, l32| within 
the picture presented here. Note that since the mean 
field theory employed here overestimates the tendency 
towards static magnetic order, our results suggest that 
the thermal conductivity in the presence of low-energy 
dynamical spin fluctuations not quite pinned by impuri- 
ties, as apparently true in YBCO, will also be suppressed 
relative to the universal limit. 

In summary, we studied the effects of disorder in cor- 
related d-wave superconductors, and calculated the DOS 
and thermal conductivity in different regimes of the Hub- 
bard repulsion U . Although the low-energy DOS is pro- 
tected due to screening of the disorder by interactions, 
disorder induced magnetism in the presence of correla- 
tions can lead to enhanced scattering which may strongly 
modify transport properties, in agreement with experi- 
ments on thermal conductivity in underdoped cuprates. 
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